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Abstract 

We present the complete electroweak one-loop corrections to top-pair production at a linear e+e^ collider in the 
continuum region. Besides weak and photonic virtual corrections, real hard bremsstrahlung with simple realistic 
kinematical cuts is included. For the bremsstrahlung we advocate a semi-analytical approach with a high numerical 
accuracy. The virtual corrections are parametrized through six independent form factors, suitable for Monte- 
Carlo implementation. Alternatively, our numerical package topf it, a stand-alone code, can be utilized for the 
calculation of both differential and integrated cross sections as well as forward-backward asymmetries. 
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1 Introduction 



At a future linear e+e collider with a centre-of-mass energy above 350 GeV, one of the most important reactions 
will be top-pair production well above the threshold (i.e. in the continuum region), 

e+ + e" ^ t + t . (1.1) 

Several hundred thousand events are expected, and the anticipated accuracy of the corresponding theoretical 
predictions should be around a few per mille. Of course, it is not only the two-fermion production process (ll.lt . 
with electroweak radiative corrections (EWRC) and QCD corrections to the final state that has to be calculated with 
high precision. Additionally the decay of the top quarks and a variety of quite different radiative corrections such as 
real photonic bremsstrahlung and other non-factorizing contributions to six-fermion production and beamstrahlung 
have to be considered. Potentially, new physics effects also have to be taken into account. For more details on 
the general subject of top physics, we refer the reader to [Ij and, for top-pair production, to the recent collider 
studies |2||3||4| and references therein. 

The electroweak one-loop corrections will be a central building block in any precision study of top-pair 
production. Also it might well be that for most of the physics a phenomenological study of the two-particle 
(top-pair) production cross section will be sufficient, thus avoiding to deal too much with many-particle final 
states observed in the detectors |5||6l0|8l|9|- For these reasons, we recalculated the complete set of electroweak 
contributions, including real hard photon corrections. Several studies on this topic are already available in the 
literature. In flO"111, the complete 0{a) corrections, including hard photon radiation, are calculated. The virtual 
and soft photon corrections both in the Standard Model and in the minimal supersymmetric Standard Model are 
determined in I12II13I . and (only) in the Standard Model in fT4l . Experience proves that so far it was difficult to 
get a satisfactory numerical comparison based on articles or computer codes without contacting the corresponding 
authors. Due to the importance of the process, for future applications, it is therefore necessary to provide a common 
basis and accessible documentation. Thus we aim, with the present write-up, to carefully document the one-loop 
radiative corrections for the process (ll.lt \ with the publicly available Fortran program topf it I17II18I . and 
with the sample Fortran outputs. In the mean time, we compared our calculations in detail with the results of two 
collaborations lil9ii20J .^ 

In this article, we sketch in short our calculation and present some typical numerical results applicable at 
typical Linear Collider energies. 

2 Conventions and Cross Sections 

In lowest order perturbation theory the process e+e^ —>■ tt can be illustrated by the two Feynman diagrams of 
Fig. 12. II For convenience we introduce the following abbreviations: 

P5=Pl+P2 = -P3 - P4, pI = t, (2.1) 

Pe = P2+P3 = -Pi - P4, pI = s, (2.2) 

P7 =P2+P4. = -Pl - P3, Pi = U- (2.3) 

In the Feynman gauge the matrix elements corresponding to Fig. 12. H are: 

'The situation witii massless fermion-pair production is much better due to the efforts related to LEP physics; see I15II16I and the 
references therein. 

Another series of numerical comparisons with the authors of [14] was started in September 2001; see also 
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Figure 2.1: Feynman diagrams for the process e+e — > in Born approximation. 



with 



M^ = — QeQt lv{p4)l^u{pi)] X [U{-P2)7f,v{-P3)], 



S-MI + iM.Tz 



Vf 



Tf-2Qf 
2 sin^w cos 



7-3 



2 sin cos 6^ ' 



(2.4) 



XP4)7'' {ve - ae75)n(pi)] X [n(-p2)7M (^t - at75)u(-p3)], (2.5) 



(2.6) 
(2.7) 



where Tj? is the quantum number corresponding to the third component of the weak isospin, eQf the electromag- 
netic charge, and 6^ the weak mixing angle. 

We parametrize the radiative corrections by means of form factors. Defining the following four matrix elements 



M\' = [v{p,)jf'G'uipi)] X [u(-p,)j^G^v{-ps)] , ij = 1,5, 
with = 1 and = 75, the Born amplitude can be written in a compact form: 

Mb = M^ + Mz = Y,. _ Fl'-^ Mt 



The form factors are 



11, B 



+ QeQt 



11, B,Z _|_ ^11,S,7 



1 ' 



^1 



15, B 



f: 



51, B 



-Ve at 



-Vt ae 



s-M^ + iM.Fz' 



s- Ml + iMzTz' 



F\ 



55, B 



Oe at 



S-MI + iM.Tz ■ 



(2.8) 
(2.9) 

(2.10) 
(2.11) 
(2.12) 
(2.13) 



Besides (12.81) . we find at one-loop level three further basic matrix element structures (in the limit of vanishing 
electron mass): 

4 

-Mlloop = E =1 5 ^a^'""' (2-14) 

a=l ''^ ' 
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with 

=-f^'G'(S)7^,&, (2.15) 

=^2G' ^ (2.16) 

Ml^ =^2G' (S) (2.17) 

= ^^^0' (2.18) 

and respectively sixteen scalar form factors Fa in total. An alternative notion uses the helicity structures, 

TWf ^ = 7^^ L 7^ R (2.19) 

etc., with L,R = (1 =p 75)/2. The interferences of these matrix elements with the Born amplitude have to be 
calculated. Only six of these interferences are independent, e.g. M^, -Ms^^ and ^ i.e. we have the 

following 10 equivalences : 

+ (2.20) 

I)Mi^^ + {s- 4m2) TWi^i _ (2.21) 
(T - [7)7^1^1 + s (2.22) 
I)Mi^^ + (s - 4m2) TWi^i - imt M3^\ (2.23) 

(2.24) 
(2.25) 

M/^^-A^4" ^ M3^^+mtMi^\ (2.26) 
Ma^^ ^ -Ma:^^ ^ M3^^ + mt Mi^^ . (2.27) 

In the massless limit {mt 0) , only Aii and Ai2 will contribute to the cross-section, and it can be expressed in 
terms of the Bom-like structures Aii exclusively. We introduce the variables 

T = ml + mf -t2:i ^{1- P cos 9), (2.28) 
U = ml + ml -uc^ cos 0), (2.29) 



47W2" 


- (T 


4 7W2^^ 


- (T 


4 7^2^^ 


- (T 




- (T 




^ 0, 




^ 0, 




^ 7W 




^ 7W 



Pt = (3 = ^l-4m2/s . (2.30) 
Based on the relations (12.201 to 12.271 ) the virtual corrections can be expressed in terms of six independent. 



^We are grateful to D. Bardin and P. Christova for drawing our attention to this simplification. 
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modified, dimensionless form factors Fl'' , F^^ , F^^ : 



^15 



?51 



?55 



F/^ + \{u-t)F^^-^iu + t + 2mi) Ff + (F4 



15 , J. o_2\ 7751 , 1/ , 7:il5 



55 



F/^ - -(u + t - 2m^) Ff + -(tx - t) F2^ 

+ ^(n - t) Ff - i(n + t + 2m?) Ff + mt (Fi^ - Ff )' 
Ff -i(n + t-2m2) F^ii + t)F 



i^s" - F, 



11 



Ff + Fi^ 



3 -^'4 

Tlie resulting cross-section formula is 
da 



p55 



p51 



mj F; 
mj Fr 



55 



15 



(2.31) 
(2.32) 
(2.33) 
(2.34) 
(2.35) 
(2.36) 



dcos 6 



™ ct P 2Ke \{u^ + + 2m?s) ( F^'F^''''* + Ff ^F= 



2s 



+ {u 



11 ^^11,-B* 
^1 

55 f.55,B* 



51 f=i51,B* 



■2 t^l 



^55^11, B* _|_ pl5pl5,B* 



+ 2mi (tn - m^) ( FgiiF,"'^* + Ff^Ff 



where the dimensionless form factors are 



f; 



,ij,B* 



s 

72 



Fl' 



mf; 



ij,B 



+ 



^ij,lloop 



(2.37) 

(2.38) 
(2.39) 



and Q = 3, a = c'^/Att. The Fa-* are defined so that double counting for the Born contributions F^*-''^ is avoided. 
The factor l/(167r^) is conventional. 

In the numerical program, helicity form factors are calculated as well. They are defined as follows: 

(2.40) 
(2.41) 
(2.42) 

i = 1 . . . 4. (2.43) 

At the end of this introductory section, we would like to give the relation of our form factors to those used in the 
literature for pair production of massless fermions. In that case, only the four form factors F^*"* contribute. They 
have to replace, in the massless limit, the form factors p and kj, which are conventionally used to renormalize 
the muon decay constant and the weak mixing angle and are precisely defined in |22 23|. We rewrite the 
matrix element A^i in such a way that it gives exactly the Bom Z amplitude (I2.5t when the four form factors 
Pet, Ke,Kt, K^t are set equal to 1: 



pLL _ 


i[ 


Fl" 


_ ^15 

I 


_ ^51 

I 


+ ^55 


pLR _ 


i[ 


Fl' 


+ Fr 


_ 7?51 

I 


-Fp 


pRL _ 


i[ 


Fl' 


_ 7?15 

I 


+ f!' 


-Fp 


pRR _ 


l[ 


Fl' 


+ Fl' 


+ f!' 


+ f!' 



E 

i,j=l,5 i,j=L,R 



Fl'M\ 
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4 e Qe at 



Pet 



(7''L 7^L) - IQel sin^ O^Ke {-/^ 7^^) 



\Qt\ sin^ (7'^L 7^) + IQt Qel sin^ 6'w«^et (t'^ 7^*) 



(2.44) 



From here, it is easy to derive the relations between the form factors F^"* in an(L,R)or (1,5) basis and pet-, Ke,Kt, Ket, 
respectively: 







4 Oe at 


S 


- M| + zM^r^ 






4 ae at 


s 


- Ml + iM^r^ 


pRL^ 




4 ae at 


s 


- M| + iM^Tz 


pRR ^ 




4 ae at 


s 


- Ml + zA/^r^ 



Pei (1 - IQel Sin^ OwKe - \Qt\ sin^ ^'w Kt + \QtQe\ Sm^ O^Ket) 

Pet {-\Qt\ sin^O^Kt + \Qt Qe\ sin^ 9^Ket) , 

Pet {-\Qe\ Sin^ 0w K-e + \Qt Qe\ sin^ OwHet) , 
Pet {\QtQe\ Sin^ dw Ket) ■ 



(2.45) 
(2.46) 
(2.47) 
(2.48) 



Three process-dependent effective weak mixing angles sin^ 9"^ and the weak coupling strength are obtained 
by inverting these relations: 



^eff = = ' " ( Fl"^ - - F^"^ + i^r) , (2.49) 

4 ae at V / 

sin^e/^'" = Ke sin^e^ = --^ — ' ^ — ^, (2.50) 



sin^ ^w'^''^'* = sin^ 6'w = —TTTIT^ ~ ' — ~ — T' (2-51) 





^pLL _ 


_ pLR _ pRL _|_ pRR^ ' 


1 




(^pRR _ pLR^ 


\Qt\ 


^pLL _ 


pLR _ pRL _|_ pRR^ ' 


1 




pRR 


\QeQt\ 


(pLL 


_ pLR _ pRL _)_ pRR^ 



(sin^C^''*)^ = Ket sii/0„ = ^ ^ (2.52) 



For the simplest approximations with factorizing, universal weak corrections, the k^, Kt, and ^/k^ become equal, 
real and constant (independent of process and kinematics); for more details see for instance the discussion of the 
weak corrections in [ 15,16 1 and references therein. There, also the important influence of higher order corrections 
is considered. 



3 Virtual Corrections 

The virtual corrections come from self-energy insertions, vertex and box diagrams, and from renormalization. 
A complete list of the contributing diagrams may be found in ifTOl . By means of the package DIANA 1241 125 1 
l26l we generated useful graphical presentations of the diagrams and the input for subsequent FORM l27ll28l 
manipulations. With the DIANA output (FORM input), we performed two independent calculations of the virtual 
form factors, both using the 't Hooft-Feynman gauge. 
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For the final numerical evaluations we used two Fortran packages: FF and LoopTools OH). Both have 
been taken from the corresponding homepages, and LoopTools was slightly adapted: one infrared Cq was added 
and the DBi was used only for photon mass A = 0. 

In the package FF, the Passarino-Veltman tensor decomposition of the amphtudes l3Tl is defined in terms of 
the external momenta of the diagrams, while in LoopTools this decomposition is performed in terms of internal 
momenta and the latter are later expressed in terms of the external ones. The resulting linear relation between the 
corresponding form factors is given in Appendix Both the ultraviolet (UV) and the infrared (IR) divergences 
are treated by dimensional regularization, introducing the dimension d = 4 — 2e and parametrizing the infinities as 
poles in e. The UV divergences have to be eliminated by renormalization on the amplitude level, while the IR ones 
can only be eliminated on the cross-section level by including the emission of soft photons. For the IR divergences 
we have alternatively introduced a finite photon mass, as is foreseen in FF, yielding a logarithmic singularity in 
this mass. Agreement to high precision was achieved for the two approaches. 

Because the calculation of one-loop corrections for two-fermion production is well known, we do not present 
a detailed prescription of the calculations. We perform the renormalization closely following |32|. On the other 
hand, we want to fix some cornerstones and sketch the renormalization and show the UV-divergent parts of all 
the contributing diagrams, such that their cancellations can be deduced. The treatment of the IR divergences will 
be discussed in more detail because of the interplay with real photonic corrections. Finally, concerning the finite 
parts, we refer to the Fortran program topf it. We only mention that we did not perform a complete reduction of 
the various scalar functions to Aq, Bq, Cq, and Dq, since this is not needed for a purely numerical evaluation. 



3.1 The self-energy diagrams 

We have to renormalize the UV singularities of the self-energies of the photon and the Z boson, and also that of 
their mixing. Since the counter terms from wave-function and parameter renormalization must exhibit Born-like 
structures, it is clear from the very beginning that a cancellation of UV divergences can only occur in terms of 
single propagator poles. The double poles, which originally occur in the self-energy diagrams, are cancelled by 
mass renormalization. This we want to stress for the following, by allowing only single poles in the self-energy 
contributions. Thus the photon self-energy and the Z self-energy diagrams take the form 

= ^ My, (3.1) 

s 

S — M| 

In the j-Z mixing diagrams, a partial fraction decomposition of the product (l/s)l/(s — M|) is performed, but 
no subtraction: 

Sz-f = -(? QtT.yz '^(jZJp ~ - ae75) X 7^. (3.4) 

We give the UV-divergent parts of the renormalized self-energies (for definitions, see Appendix iBt: 



23\ 1 

o2 



= = 2 Sing. COS U « - - Y') 7' ^'-'^ 
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V^ren.UV 2 / s^2n / , ^ 41 \ 1 

= e (g - M^) — + . - _— -. (3.7) 

V o sm 6 cos"^ 6^ sm t^^ / ^ 

Three families of fermions are assumed. The UV-divergent terms of the self-energies are independent of the 
fermion masses, but this is, of course, not true for the finite contributions. 

The form factors are easily deduced from the above representations. The photon self-energy, for instance, 
contributes to F^^ only: 

yren 2 

Pll,77 ^ ^11,5,77 ^ = Q^Q^^ (3.8) 



3.2 The vertex diagrams 

From the initial-state vertex corrections, form factors -F^'^'^, V = j,Z, arise, and from the final vertices F^-''^ and 
F^-''^, the latter being proportional to m^. There are UV divergences from the vertex diagrams in F^-': again only 
Bom-like amplitudes are UV-divergent. 

The divergent parts of vertices with a photon or Z boson in the s-channel, correspondingly, are:^ 



7 



Y 1 1 



2sw 



S € 
2 



fi'^^it. ^1"+ ® 1^15 + /r'^7M75 
1 1 



pl5,7. 



f51.7^ 



Mi e 



The expUcit expressions from the initial and final photonic vertices are: 



(3.9) 



'7''75 



(3.10) 



11,7 _ 



n 



17 1 



64 
27 



1 mi 



5 1 

-3 



2 

+ 3 



32 
¥ 
10 

y 



5 1 



For the Z boson in the s-channel we only give the sum of the initial- and final-state fermion vertices: 



(3.11) 
(3.12) 
(3.13) 



11,Z 



/i 



15, Z 



fl 



fl 



55,Z 



973 
216 

21 



665 
216 
97 
72 " 



25 1 
18^ 



9 rrii 



1 

16 m|; 



2 "T-t 



1 2 ]n 

3^ Ml 



1 



157 2 82 4 3 



7 ^ 1 



-w 



16 Ml 



16 Ml 



■ 2 "f^t 



1 



2 

Ml + 3^^^" 



108 
137 



z 

(3.14) 



32 



2 ^ 5 T?X^ 



1 ml 



3 Mr 



95 1 
"108 4; 

7 1 
1272" 



16 M2, 



7 mf 



+ 



1 ml 
~ 48M^ 
1 ml 



449 



10 



1 



1 mr 



lOs""^^ 3"^^ AMl^ 12 Ml' 



77 



1 m} 



48 M2 48 M2 



36"''^ ^ 12 M2 



1 ml 

umI' 



The resulting form factors can be extracted, e.g. 



»J>7 



(3.15) 
(3.16) 

(3.17) 
(3.18) 



''To compactify the following formula we introduce the abbreviations: sin 9^, = sw and cos Oy, = cw. 
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3.3 The box diagrams 



ZZ, Z7 and 7Z box diagrams contribute to all form factors F^^ to F'^ introduced in (I2.14t to (I2.18t . while the 
WW box diagram contributes only to F^^ . The pure photonic box diagrams contribute only to F^ and F^ . 
Simple power counting shows that there are no UV divergences from the boxes. The IR divergences will be 
discussed in Section l331 and AppendixO 



3.4 The counter- term contributions 

Finally we have to take into account the contribution from the counter terms of Appendix |b1 where we also 
introduce some of the notation to be used in the following. With these the photon exchange becomes: 



be 



QeQt — • 
S 



(3.19) 



Analogously, for the Z exchange 



Cz 



^lj.{Ve - aejb) ^ liJ.{vt - atj5){Za,t + Zb,tlb) + liiiVe - ae75)(2a,e + Zb,el5) ^ Itiivt - atj^) 



7m ® 7m {vt - 0*75) 



sw cw 



OSiv - Ji,{Ve - ae75) O 7;^ OS 



Sw Cw 



w 



(I 1 w 2 



s-Mf 



(3.20) 



It is again easy to collect from the above expressions the corresponding contributions to the form factors F^^ to 
Ff^. The contributions to, say, F^^ from the counter terms are: 



^ll,ct 



Je 

2 h Za.t + Za, 

e 



{vtQe + VeQt) 



OSw 



e 



2— + Za,t + Za,e + —] ^Sy^ 



1 



1 



F, 



ll,B,Z 
1 



Sw Cw s — M| 



The resulting 1/e terms may be read off in Appendix IbI 



(3.21) 



The sum of all the 1/e terms listed in the foregoing subsections for the form factors Fl\ i,j = 1, . . . , 4, has 
been shown to finally vanish separately for the photon and the Z pole of the s-channel propagator: 



pij,VV _ pvm _|_ _|_ pij,Zj _|_ pij,ZZ _|_ pij,j _|_ pij,Z _|_ pi 



ij,ct 



UV 



0. 



(3.22) 



3.5 Infrared divergences 



In topf it, we have incorporated two weak libraries. One uses the package LoopTools l30l . and the other one 
the package FF |29|. With these two options, we have a variety of internal cross checks at our disposal. 

The photonic virtual corrections contain infrared divergences. They appear as singular behaviour of classes of 
one-loop functions. One may follow several strategies to handle them in a numerical calculation. The simplest one 
is to blindly give the task to the library for numerical calculation of the one-loop functions and then control the 
infrared stability numerically in the Fortran program. Both packages allow for this approach; LoopTools with 
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dimensional regularization or with a finite photon mass, while FF treats loop functions with finite photon mass 
only. 

In addition, we checked the IR stability in two ways explicitly. In the library based on FF, we simply took a 
small but finite photon mass and directly applied FF without simplifying any tensor functions. Several analytic 
checks were also performed. In the other one, we isolated in all the scalar functions the IR divergence explicitly 
and the cancellations with the divergences from bremsstrahlung corrections (see Section l43l were controlled both 
analytically and numerically. 

In AppendixO we fix the notation and discuss the basics of the treatment of IR divergences. In the rest of this 
section, we simply give a list of the divergent parts of the form factors: 

From the renormalization of the fermion self-energies (wave function renormalization factors) : 

pij,Zf,m ^ -4e^Q}m} DBi{mf,m},0) Fi^^"". (3.23) 
From the vertex corrections (index / = e, t): 

pij,Vf,m ^ -2e'^Qj{s-2mj)Co{m},s,mj;0,mj,m}) Fi^^^. (3.24) 

These form factors combine in the cross-section with the initial- and final-state soft photon corrections to an 
infrared-finite contribution. For instance, the pure photonic parts contribute only to F^^. The resulting IR-divergent 
cross-section contribution in (I2.37t . 



da^'^^ da 



B 



a 



2 xlR 



(3.25) 



6f = 2ln^ I 1 + ^ In | , (3.26) 



d cos 9 d cos 6 n 

with 



.IR _ o 1. 'IV 

A 




is compensated with (14.761) and (I4.77t . 

A little more involved are the box diagram contributions. As a typical example, we show the photonic box 
parts. The direct box gives: 

pii,dj,m ^ _p55,di,m ^ ^t_^_s^g^^ (327) 

^n,d7,iR ^ ^55,d7,iR ^ _4g^^ (3 28) 

^55,d7,iR ^ -AmtQd, (3.29) 



with 



Qd = -e^ QeQt Co{mlt,ml,0,mlm^) fI^''''\ (3.30) 



From these expressions, the form factors (I2.31t to ( I2.35t get contributions: 

pii,d7,iR ^ +AT QeQt Co{mlt,ml,0,mlmf) fI^'^'^, (3.31) 
and all the others vanish. Analogously, from the crossed photonic box diagram: 

^ll,C7,IR ^ ^^155,C7,IR ^ _2(^_^2)^^^ (3 32) 

^n,C7,IR ^ _^55,C7,IR ^ _^g^^ (3 33) 

^55,c7,iR ^ +4mtgc, (3.34) 
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with 



-e^ QeQt Co{ml,u,mf -,0,7711,1711) fI^'^'"' . (3.35) 



From these expressions, the form factors (I2.31t to (I2.35t get contributions: 

^1M7,IR ^ -AU QeQt Co{7Jil,u,7Jil,Q,ml,ml) fI^'^''^ , (3.36) 

and again all the others vanish. 

The resulting cross-section contributions become 

4 - QeQt In- In^^ - In^-^ , (3.37) 



d cos 6 d cos 9 TT X \ T U 



and are compensated with (I4.78t . 

In Appendix Owe show the relevant scalar functions explicitly. 

4 Real Photonic Radiative Corrections 

4.1 The three-particle phase space 

The reaction 

e+(p4) + e-(pi) ^tfe) + tfe) + lip) (4.1) 

with 

is the one which in reality always takes place, even if for soft photons the 'elastic' Born cross-section can be a 
good approximation without taking into account the radiated photons. Here we introduce the final-state phase- 
space parametrization: for convenience of notation, the top physical momenta q2 = —p2 , = —ps are used here. 
We will not neglect the electron mass systematically, p\= p\ = mg, and (5o = y^l — Ami/ s. ^™ semi-analytical 
integration approach with physically accessible observables as integration variables may be used to set benchmarks 
and to control the numerical precision to more than four digits. Their choice is constrained by the observables we 
want to predict, notably the angular distribution and certain cross-section asymmetries. Basically we follow the 
approach proposed in 13311341 and extend the required formulae to the massive fermion case. 

There is not too much found in the literature for the radiative production of massive fermion pairs. Thus, we 
will discuss the kinematical details with some care, since they define the integration boundaries of our numerical 
integration program. 

In (12.11) and (12.31) we defined t and u for two-particle production. With the additional photon in the final state, 
we have to be more specific and will use the following definitions: 

T = 2pm, (4.3) 
U = 2piq3. (4.4) 

Additionally, the following invariants will be used: 

s' = (92 + 93)', (4.5) 
Zi,2 = 2ppi,4, (4.6) 



11 



Vl,2 = 2pg2,3- 



(4.7) 



The squares of all three-momenta in the centre-of-mass system can be expressed in terms of a set of A functions: 

4s|g2|2 = Ai = \[{pi+PA)\{q^+P?,ql] = {s-V2f -Amis, (4.8) 

Asm'' = ^2 = X[{pi+P4)Mq2+pf,ql] = {s' + V2f-4mls, (4.9) 

As\pif = 4s\p4f = Xs = \[{pi+P4f,plpl] = s'-4mls, (4.10) 

As\p\' = Ap = X[{pi+PA)Mq2 + q3f,p''] = {s-s'f, (4.11) 

where we use A (x, y, z) = x^+y^+z^ — 2(xt/+a;z+yz) and the relation 4p^|p*Bp = X[{pa+Pb)'' ,p\,p'b]\pa=o 
for PA = {pi + Pa) = {^/s, 0, 0, 0). 

The phase space of three particles in the final state is five-dimensional. This means that only five of the ten 
scalar products (those introduced already: s, T, U, s' , Zi, Vi, plus Wi^2 = '^Pi,Aq2) built from the five momenta are 
independent. In fact, the following relations hold 

s = s' + Vi + V2 = s' + Zi + Z2 = V2 + Wi + W2 (4.12) 
= Vi + U + T = Zi + Wi + U. (4.13) 



We use the first of relations (I4.12t in order to substitute Vi everywhere in favour of V2 as already done in 
(I4.9t . Since we do not consider transversally polarized initial particles, an integration of the cross section over 
the corresponding rotation angle is trivially giving a factor 27r and we are left with four non-trivial phase space 
variables. 

For the calculation of the forward-backward asymmetry, the angle 6 between the three-momenta of t and e+ 
is used, in accordance with (I4.3t . Further, the energies Et, Ei, and E^ are 'good' variables. As mentioned before, 
they can be expressed in terms of the invariants s' and V2: 

E-y = (4.15) 



2^5 

Et = ^rs-E^-E-t = ^-^^. (4.16) 



V2 
2^s ■ 



The scattering angle in the centre-of-mass system may now be expressed by invariants: 



cose. (4.17) 

As will be seen later, the two invariants s' and V2 also describe the angles between any pair of final-state particles. 
Therefore, we choose them to parametrize the phase space. Finally, the fourth integration variable will be the 
azimuthal angle of the photon 

The coordinate system is chosen such that the i moves along the z axis and the beam axis is in the y-z plane. 
The four-momenta of all particles can then be written as follows: 

pi = ^(l,0,-/3o sine, -/?o cos 0), (4.18) 

PA = ^(l,O,/3osin0,/3ocos0), (4.19) 
p = Ej{l,sm9^cos(j).y,sm9^sm(j)^,cos9^), (4.20) 
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qs = {Et,0,0,\qs\), (4.21) 

92 = Pi+Pi-p- 93- (4.22) 

The (p'y and 9^ are the azimuthal and polar angles of the photon. The expression for cos 6^ (and also that for jgsl) 
can be obtained from (p + 9*2)^ = (?3)^> 

cos^, = ^--^^-^'''''7''-^^+''', (423) 
and again depends only on s' and V2. The differential bremsstrahlung cross section (I4.2t takes the form 

^ • — d(^^drda;dcos6'. (4.24) 



(27r)5 2s/3o 16 
In the last step, dimensionless variables are introduced: 

V2 

X = —, 
s 

s' 

r = -, 
s 

1"m, — ■ 
S 



(4.25) 
(4.26) 
(4.27) 



The integration boundaries are either trivial (0^ and cos 6) or can be found from the condition that the three three- 
vectors p, q2, 9*3 form a triangle with the geometrical constraint cos^ 6^ < 1. The four integration variables vary 
within the following regions: 

< (p^ < 2tt, (4.28) 

X 



1 + x- ^y{l-x)'^ -rm) < l-r <— 1 + X + V(l - - r„ , (4.29) 



2x + rm/2 \ V V . my - - 2x + rm/2 

< X < 1 - (4.30) 
-1 < cos 6* < +1. (4.31) 

If the order of integrations over r and x is interchanged, their boundaries are 



'-^^-Jl-^]< X <i^fl + Jl-^), (4.32) 



r J 2 \ \ r 

rm< r < 1. (4.33) 
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The kinematic regions of r and x are shown in Fig. l4.ir a) for massless (r^ = 0) and in Fig. l4.ir b) for massive 
{Tm / 0) final fermions. At the kinematic boundaries, the three-momenta p, 5*2, 9*3 are parallel. Further, there are 
three special points where exactly one of the three three-momenta vanishes 



A 

B 
C 



(1,0). 



1 



:, (1 - Vr^) 



(4.34) 

(4.35) 
(4.36) 



At point C the soft photons are located. Section |431 is devoted to their treatment. The t (t) are at rest in A (B). 
In the massless case, = 0, the three points A, B and C are located at the comers of the kinematic triangle, 

(1,0). 



A = (0, 1), B = (0, 0), C = (1, 0). From (I4.14t - (l4.16t it follows that the photon energy is maximal at the left 
edge, coinciding with the x-axis; the fermion energy is maximal at the lower edge, coinciding with the r-axis; and 
finally the energy of the anti fermion is maximal at the third edge. 



4.1.1 Energy cuts 

Cuts on the energy of the final state particles are of importance for two reasons: they are being applied in the 
experimental set-ups, and for the photon we have to identify the soft photon terms in order to combine them with 
virtual corrections for a finite net elastic cross section. The lower hard photon energy (being also the upper soft 
photon energy) is 

io = E^'''. (4.37) 
All three energy cuts are deduced from (I4.14t to (14. 16t . The photon energy is related to r by 

r = 1 - 2E^l^s (4.38) 
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Figure 4.2: Triangle of the three-momenta of fermion, anti fermion and photon. 



and the limits to be imposed are 

r^ = l- 2E°^^^/Vi < r < 1 - 2uj/^ = r^. (4.39) 
Constraining the fermion energies leads to cuts on x: 

r^-r = 2Ef''/^fs -r <x<l-rE = l- 2Ef'''/^/^. (4.40) 

All the energy cuts are independent of the mass of the final fermions and of cos 9. They are illustrated in Fig |4.1l 
From (14.401) it can be seen that the derivatives of the kinematic border at points A and B in Fig. I4.1f b) (for their 
definitions see (I4.34t and (I4.35t ) are and —1. 



4.1.2 Angular cuts 

The scattering angle 9 

is the angle between t and . This angle is one of the integration variables and is constrained directly: 

Cmin < C = 008 6* < Cmax- (4.41) 

Additional angular cuts deserve a study of the (r, x) parameter space. To be definite, we will always consider the 
kinematic bound of r for an arbitrarily chosen value of x. 

The directions of the final-state particles define three angles 9ii, 9t'y, and 9^^; as shown In Fig. 14.21 The 
acollinearity angle ^ is defined by 

e = TT - 9u. (4.42) 

The condition ^ ^ 1 restricts the events to a Bom-like kinematics: the fermions are back to back and only soft 
photons or photons coUinear to one of the final fermions are allowed. Using the above formulae, the acollinearity 
angle ^ is expressed in terms of the invariants x and r: 

cos C = ^^+;^^^ = r(l + x)-x(l-x)-r^ _ 

2\/XiX2 V (1 — x)2(r + x)2 — rm[(l — x)^ + (r + x)^] + 



' m 



Equation (14.431) can readily be derived considering the scalar product 52^3 or alternatively the triangle of the three- 
momenta of the final particles in Fig. 14.21 with account of the relations (14.111) between A functions and absolute 
values of three-momenta. 
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0.2 0.4 0.6 0.8 ^ 1 0.2 0.4 0.6 0.8 ^ 1 



Figure 4.3: The kinematic region of r and x for different values of the acoUinearity angle ^ for (a) = 4mf /s = 
and (b) = 0.2. 



For massless fermions, Eq. (I4.43t is much simpler and describes a hyperbola with a symmetry axis rotated by 
an angle of — vr/S relative to the r axis. The kinematic regions for different values of the acoUinearity angle are 
shown in Fig. l4.3l All lines intersect at the points A and B. For moderate cuts on the maximum acoUinearity angle, 
the kinematic area is only constrained for values of x above the point B, i.e. for x > y^rv^(l — y^r^) /(2 — ^/r^). 
In this case only the lower bound of r is changed. The constraint to the kinematic region acts in a way similar 
to a cut to hard photons. This is clear from the topology of events with high acoUinearity: the fermion and anti 
fermion fly approximately in one direction and must recoil against a hard photon. For stronger acoUinearity cuts, 
constraints of the kinematic area arise also for values of x below the point B. In this case the allowed range for 
r splits in two regions. The first region extends from the lower kinematic border to the smaller solution of (I4.43t . 
while the second region extends from the larger solution of (14.431) to the upper kinematic border. 

The analytic treatment of the acoUinearity cut for the massless case was also discussed in l35U36ll23ll37l . 

In a similar way as explained for the acoUinearity angle, the two angles 9tj and 9^^ can be expressed in terms 
of the two invariants r and x: 

cosO,, = ^^-^"^= r(l + .)-(l-x) 

^ 2y^p (l-r)V(l-x)2-r„' 

n - n Ai - Aa - Xp -x(l + r)+r(l-r) 

cos6t^ = cos 6^ = , = , (4.45) 

^ 2v/A^ (l-r)V(x + r)2-r„ 

Although physically the situation for a cut on 9tj is equivalent to a cut on 9[^, the symmetry is broken because we 
had to make a choice between Vi and Va- The constraint on 9tj leads to a quadratic equation in x and r, while 
that for 9i^ is quadratic in x and of fourth order in r. In the massless case, the constraints to 9fy become a bilinear 
equation in r and x, describing a hyperbola with a symmetry axis, which is rotated by an angle 7r/4 relative to the 
r axis. The massless limit of ( 14.451 ) leads to a constraint on 9i^, which is linear in x and quadratic in r, describing 
a hyperbola with a symmetry axis, which is rotated by an angle of 7r/8 relative to the x axis. 

The kinematic regions for different values of the angle 9t'y are shown in Fig. 14.41 All lines intersect at the points 
A and C. The exclusion of events with small angles 6f^ excludes regions near the edge x = 1 — r. These kinematic 
regions correspond to events with large anti fermion energies, (compare Fig. 14. II) . Technically, a constraint of the 
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angle Ot^y from below changes the upper bound of r in the kinematic region for a fixed x. The lower bound of r is 
unchanged by this cut. 

Finally, the kinematic regions for different values of the angle 9^^ are shown in Fig. 14.51 All lines intersect 
at the points B and C. The exclusion of events with small angles Oj^ excludes regions near the r axis. These 
kinematic regions correspond to events with large fermion energies, (compare Fig. l4.2t . Technically, a constraint 
of the angle 6^^ from below affects the kinematic region only for some x below a certain value. For these x and 
for massless fermions, the integration region of r is split into two regions. The first region extends from the lower 
kinematic border to the smaller solution of (I4.45t . and the second one extends from the larger solution of (I4.45t 
to the upper kinematic border. For massive fermions, the cutting out of small angles %^ changes only the lower 
bound of r as far as x is below the point B. For a harder cut on 9i^, the kinematic region is also affected for values 
of X, larger than the x coordinate of the point B. For a fixed x and finite mt, two regions of r near the kinematic 
border are then allowed, while some region in the "middle" is cut out by the constraint. This is similar to the 
situation in the massless case. 



4.2 Radiative differential cross sections 

For massless fermions, typically a threefold analytical integration of the radiative contributions to fermion pair 
production with realistic cuts may be performed, see 136, 23l l35ll and references quoted therein. For massive 
pair production, everything becomes non-trivial and, in the end, we decided to perform only the first integration 
analytically, that over 0^. This leaves three integrations at most for a numerical treatment. Our practice proved 
that the accuracy and speed are absolutely satisfactory for our needs of calculating benchmarks and physics case 
studies. 

For this reasons, and in order to make everything well-defined, we now have to collect the singly analytically 
integrated contributions to be used in a subsequent numerical calculation. 

We will use a notation for the couplings with some flexibility not needed in the Born case: 

V{S,S') = \^er^[xMX*n{s')+Xrn{s')x*n{s)] 
\ m,n 
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Figure 4.5: The kinematic region of r and x for different values of the angle and for (a) = 4m| /s = and 
(b) rm = 0.2. 



X [ve{m)vl{n) +ae(m)a*(n)] • [vf{m)vf{n) + af{m)a*f{n)] j, 

Ais,s') = l^elj2i^nr{s)xUs')+Xrn{s')xUs)] 
(, m,n 

X [fe(m)a*(n) + ae{m)vl{n)\ ■ [vf{m)a*j{n) + af{rn)v*j{n)\^, 

C{S,S') = 2^e}^[Xnr{s)x*n{s)+Xnr{s)x*n{s)] 
\ m,n 

X [ve{m)vl{n) + ae{m)a*^{n)\ ■ af{m)a*j:{n)^, 

C{S,S') = ^el^[Xrn{s)x*nis')-Xm{s')Xnis)] 
[, m,n 

X [ve{m)al{n) + ae(m)f*(n)] • [vf{m)a*j:{n) — af{m)v'j^{n)]^, 
where we use Vf{^) = Qf, 0/(7) = 0, Vf{Z) = vj, and af{Z) = a/ and 

^^^'^ " s-Ml + iFzMz' 

With these conventions, the Bom cross section becomes 



(4.46) 



dc 



2s 



1 - /?2 

V{s, s){2 -0^ + c^0^) + 2c(5A{s, s) s) 



(4.47) 



(4.48) 



(4.49) 

(4.50) 
(4.51) 

(4.52) 
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The cross section for e+e" — > tfy subdivides in the gauge-invariant subsets of initial-state radiation, final-state 
radiation and the interference between them. Explicit expressions for the totally differential cross section may be 
found in 137.1 . The integration over 0^ is not too complicated and in fact we could simply use existing tables of 
integrals I33[l34ll37ll . This first integration is unaffected by the cuts discussed and has to be performed with an 
exact treatment of both rrig and nit. 

The cross section for initial-state radiation after the integration over c/)^ is 



dcds'dV2 



with 



238'"^ 

-2 



V{s\s') 



2^li%(2r2 - 2Ts' + s'2 + 2mis') 



„'2 



+ 



VD2 

2s' 



S — S' 

1 



,^{2U^ -2Us + s'^ + 2mis') 

= + -L^ - Ts' + C/2 - Us' + s'2 + 2mls') 



1 



{-2Us' + ss' + + 2ml{s + s')) + —={-2Ts' + ss' + s'^ + 2ml{s + s')) 
Di V-C2 



-2s' - 4mt 



+A{s',s')s 
2s' 



s — s' 



1 1 

+ 



^07 



+C{s',s')mf 
-{s + s') 



2s'ml 



(U-T) 



Ci I C2 



1 



1 



(s + s' - 2U) + ^=(s + s' - 2r) 



2s' 



/2 



1 1 

+ 



+ 2 



s — s' 



+ 



1^2 
1 



(4.53) 



Cl,2 

/3o 



Dl^2 — ■T-Cl,2 
A2 



s — s 



(1 ±Pocct) , 



yi(s + s')-s(s-s') 



(s - s')V>^ 



4m^ 



1 f 1 r 



+ Ami 



2ss' - {V2 + s')(s + s') ± c/3o(s - s')VA2 
s'y2(s - s' - F2) -{s- s'fm^'j } . 



(4.54) 
(4.55) 

(4.56) 
(4.57) 



The cross section for final-state radiation after the integration over (pj 

d^- 



IS 



dcds'dV2 



2s3 



Vis,s) 



-'^^{-2UT + s^ + 2mis) 
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2m| 



_^(_2(C/r + UC2 + TCi + C12) + s' + 2mis) 



( U{T + C2) + T{U + Ci) - 2mf s) - |^(2C/r + 1(2 + TCi + C12 - 25^) 



+ (^2 _ ^rni) + —{V,-2s- 4mf) 

Vi V2 



+A{s,s)s 
1 



HiT-U) + ^-^{T + C2-U-Ci) 



f^(.'-2m?)(2C/ + Ci-2r-C2) + ^^ + ^ + ^' ^ 



V^1^2 



^2 



+C{s, s)2ml mis (^1^ + 1^+ :^(Ci2 - ss' + 2mls) | , 



(4.58) 



with 



C12 = C1C2 - i(s - s'fplil - c'){l - cf). 



(4.59) 



The cross section for the interference between initial- and final-state radiation after the integration over is 



dcds'dV2 



2s2s 

1 



V{s,s') 



{s - U){2T^ - 2Ts' + s'^ + 2mis' - 2UT + - 2m|s) 



+ 



D2V1 
1 



-D2F2 
1 



{s - T){2U^ - 2Us' + + 2mls' - 2UT + - 2mls) 
{s' - U){2T^ - 2Ts + + 2ml s - ^UT + s'^ - 2ml s') 



(s' - T){2U^ - 2Us + + 2mj2s - 2UT + s'^ - 2mf s'^ 



+pr ((Ci - C2)s + {U- T){3s -s'- Ami)) 

((Ci - C2)(^' + 4m|) + {U- T){3s' - s + Ami)) 



+ 



1 



+A{s,s')s 
1 



D2 VDl 
1 



s^ + s'^ + 2mt2(s + s')) 



{s - U)[-2Ts' + s'^ + s(C/ - T)] 



+ 



D2V1 
1 



(s - T)[-2Us' + s'^ + s{T - U)] + 
{s' - T)[-2Us + s^ + s'{T - U)] 



D2V2 



[s' - U)[-2Ts + s^ + s'{U - T)\ 



{2ss' + 2ml{V2 - 2s)) + i-2ss' + 2ml {Vi - 2s)) 
Vi V2 



1 



1 



+^={s'T - sU) + -^{s'U - sT) + 2{s + s') + Ami 



+Cis,s')ml (s + s') 



^ {s-U) + ^^{s-T)-^^rTT- 



D2V1 



D2V2 



D1V2 



-U 
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' Vi 

+C{s, s')'mf{s + s' 



1 1 



(4.60) 



4.3 Soft photon corrections 

The four-dimensional integration of the bremsstrahlung contributions is divergent in the soft-photon part of the 
phase-space and is treated in d dimensions. One starts from a reparametrization of the photonic phase-space part 
with Born-like kinematics for the matrix element squared. To obtain a soft photon contribution we have to take the 
terms of the bremsstrahlung amplitude without = Ej < lu in the numerators. In this limit, s' approaches s and 
the soft contribution to the differential cross section takes the form 



decs 6 



a 



J Born 
_§soft ^ 

vr d cos 9 



(4.61) 



with 



47r^ 
1 

47r 



d^p 



{2TTf2E^ 
d^p 



7 



Qe 



2P4 
Z2 



2pi 



+ Qt 



'202 
Vi 



203 

V2 



e{uj - E^) 



and 



jsoft 

4^2" 



E^)r°f* 



2mt 



(4.62) 



Z\Z2 

s 



+ QeQt 



+ 



u 



u 



ZiVi Z2V2 Z1V2 Z2V1 



2m^ 



^2 



V1V2 



(4.63) 



The scalar products have to be taken according to Bom kinematics, i.e. the expressions (I4.6t and (14.71) become 

Zi = 2ppi = 2E^ [pI + \pi I cos 9p\ , (4.64) 



'Pi ' 

Z2 = 2pp4 = 2E^ [pI - \pi\ cos 9p] , 
Vi = 2po2 = 2E^ [gO + \o2\ cos 9q] , 
V2 = 2po3 = 2E^[oi-\o3\cos9q]. 



(4.65) 
(4.66) 
(4.67) 



From here we see that /*°-^* is constructed to be independent of Ej. Substitute now, with d = 4 — 2e, e < 0, 



^soft 



(2vrM) 
47r 
1 rl 

2 



2e 



E~^^+^^)dE, 



di 



UJ 



1 



TSOft 



We introduce the abbreviation for the infrared divergence 

VlR 

The infrared divergence can also be regularized by introducing a finite photon mass A: 

PiR - In/i 



+ f -l„(2Vi). 



'4- 



(4.68) 



(4.69) 



(4.70) 



21 



The last integral over ^ = cos 9^ is trivial for the products ZiZj and ViVj, since they contain only one angle; 
one may thus identify either ^ = cos Op or ^ = cosOq. In the initial-final interference, one may introduce a 
Feynman parameter 



1 



1 1 



1 



ZiVj 4 ppi pqj 4 



da 



1 



with 



(4.71) 



kij = api + (1 - a)qj. 
Then, (pkij)'^ = E^s{l — [5ij cos 6ij)'^/4: and identify now ^ = cos 9ij. Further, 



(4.72) 



(l - Puf = {I - f^Asf = {I - Prf = -[ail-a)T + a'^ml + {l-afmf], (4.73) 

s 

(1 - Pnf = (1 - p^2f = (1 - Pu? = - [a(l - a)U + a^ml + (1 - afn?^ . (4.74) 



The final result is 



with 



and 



2 [PiR + ln^-^) 



s - 2m; , 
-1 H ^-^In 



s/3o 



l-/3o 
+ Li2(l) +Li2 

xsoft 



s(3o 



l + /?o 
l-/3o 

2/?o 
l-/3o 



/3o 



2/3o 



+ Li 



2/3o 



(mt,cj,e,/x), 
2(P,^ + ln^l 



In 



/5o + l 



in ^= 



1 



+ -[TJ^{T)-UJ^{U)], 



da- 



In 



1 + 

1-/37 



(4.75) 



(4.76) 
(4.77) 

(4.78) 

(4.79) 
(4.80) 



and analogue definitions for T ^ U. We calculate the finite interference part given in ( I4.78t numerically, but have 
shown the agreement with Eq. (3.64) of l38ll : 



T T{T) - U T{U) 



Lio 1 



Lis 1 



1-P 



1 - ^ cos 61 

l-(3 
1 + /? cos 6* 



+ Li2 1 



Uo 1 



1 + 



1 - PcosO 
I + (5 cos 9 



(4.81) 
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5 Results 



In this section we present the numerical results of the electroweak one-loop calculation to the process e+e" tt . 
We have performed a fixed-order a calculation, i.e. no higher order corrections such as photon exponentiation 
have been taken into account. 



For the numerical evaluation we assume the following input values I20lll9irr8l : 

= 2.49977 GeV , a = || = 1/137.03599976 , S:^^^ = ^/10^ 

Myy = 80.4514958 GeV , Mz = 91.1867 GeV , = 120 GeV , 

me = 0.00051099907 GeV , mt = 173.8 GeV , mb = 4.7GeV , (5-1) 

m^, = 0.105658389 GeV , m„ = 0.062 GeV, md = 0.083 GeV, 

rrir = 1.77705 GeV , rUc = 1.5 GeV , m, = 0.215 GeV . 

Two packages, namely FF f29l| and LoopTools |30| have been used for the numerical evaluation of the loop 
integrals. 

In Fig. 15. II we present the differential cross section for various generic values of ^/s. 



1 





Figure 5.1: Top-pair production: Differential cross sections in Bom approximation (solid lines), with full 
electroweak corrections (dashed lines), with an s' = 0.7 s cut (dash-dotted lines); also shown: pure weak 
corrections (dotted lines, photonic corrections and running of a excluded); all for ^/s = 0.5, 1, 3 TeV. 
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It can be seen that for rather high centre-of-mass energies the characteristic features of a massive fermion 
pair production become less prominent. At ^/s = 3 TeV the differential cross section of electroweak radiative 
corrections starts to exhibit collinear mass singularities at the edges of phase space. Those are cured by applying a 
cut on s'. In general it can be seen that the effects of radiative corrections are more dramatic for top-pairs produced 
close to the direction of the beam. For the TESLA range of centre-of-mass energies, backward scattered top quarks 
give rise to slightly larger corrections than forward scattered ones fl9]. For higher energies this effect is more or 
less washed out. 

In Table l5.ll to 15. 31 we present a complete set of form factors entering the cross-section calculation. The form 
factors given correspond to the minimal set of independent form factors possible for a two-to-two process with two 
massless and two massive fermions in the initial and final state respectively, and are defined with respect to the 
'naturally' arising form factors in Eq. (I2.31t . For completeness we also give the corresponding Bom form factors. 
The numerical values given are obtained for a characteristic centre-of-mass energy of ^/s = 500 GeV and a fixed 
scattering angle cos 9 = 0.7. 



f.f. 


Born 

Re Im 


weak 1-loop contributions 
Re Im 


p51 
Ff 

mt Fg" 
mt 


-2.5092647 lO"'^ 6.0265891 lO^^^ 
1.562008310"*^ -1.473211910"^^ 
5.623996310"® -5.304285710"" 

-1.374797210"^ 1.2966433 10"^° 
0.0 0.0 
0.0 0.0 


1.1805990 10"® -3.2896119 10"9 
-1.050791510"® -8.462730310"^ 
-7.7050611 10"^ -5.8986660 10"^ 
-4.8821798 10"^° 7.4750196510"^ 

9.088270510"^° -8.9067902 10"^° 
-9.5315102 10"^° 5.099511710"^° 



Table 5.1: Real and imaginary parts of the six independent form factors F/*^ = ^ F^^ for weak, non-photonic 
corrections to the process e+e" —>■ ti at ^/s = 500 GeV for a fixed scattering angle cos 9 = 0.7. For reference we 
also give the corresponding Born form factors. 



f.f. 


Born 

Re Im 


weak 1-loop contributions 
Re Im 


pl5 

p55 

mt F3" 
mt Fg^ 


-6.269143510"® 3.5795119 10"^^ 
3.806796210"^ -8.7501891 lO"^^ 
1.370633410"® -3.1504975 10"^^ 

-3.3505411 10"® 7.7014546 10"^^ 
0.0 0.0 
0.0 0.0 


5.189045910"° -3.119828110"^° 
-3.790426010"° -2.418045110"° 
-3.025350910"° -1.856707610"° 
1.0657621 10"° 2.3331039 10"° 
1.2278260 10"i° -5.840895310"" 
-9.930667210"" 4.245246910"" 



Table 5.2: Same as TableOfor ^/s = 1 TeV. 
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f.f. 


Born 

Re Im 


weak 1-loop contributions 
Re Im 


pl5 

p55 

mt F^^ 
rrit Fg^ 


-6.9644350 10"^ 4.3540072 10"^^ 
4.1984821 10-1° -1.0643459 10"^^ 
1.511659610^9 -3.8321674 10"^^ 

-3.695282310^9 9.367810410-1-^ 
0.0 0.0 
0.0 0.0 


1.015016210^9 5.64937310-12 
-6.752650810-1° -3.423640010^1° 
-6.075180810-1° -2.675495810-1° 
3.563240010-1° 3.497406710-1° 
2.989516310-12 -6.6708986 lO'i^ 
-2.493916010-12 9.1292861 lO-i^ 



Table 5.3: Same as Table 15.1 I for ^/s = 3 TeV. 

In Fig. l5.2l we present the total integrated cross section as a function of ^/s. From the previous discussion it is 
clear that the effect of radiative corrections is less dramatic in the total cross section, since the effects above and 
below the Bom cross section are averaged out. 
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Figure 5.2: Total cross-section for top-pair production as a function of s. Born (solid lines), electroweak (dashed 
lines), electroweak with s' = 0.7s-cut (dotted lines) and electroweak with s' = 0.7 s- and cos 9 = 0.95-cut 
(dash-dotted lines). 

Finally the forward-backward asymmetry of the total integrated cross section can serve as a good observable 
to determine the effects of radiative corrections. Towards higher energies, the effects become distinctively. 
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0.7 




Figure 5.3: Forward-backward asymmetry for top-pair production as a function of s. Bom (solid lines), 
electroweak (dashed lines), electroweak with s'-cut (dotted hnes) and electroweak with s' = 0.7 s- and cos 9 = 
0.95-cut (dash-dotted hnes). 



In summary our calculation shows that for the next generation of linear colliders with centre-of-mass energies 
above -y/s = 500 GeV, electroweak radiative corrections modify the differential as well as the integrated cross 

section within the experimental precision of a few per mille. The package topf it provides the means to calculate 
those corrections and allows predictions for various realistic cuts on the scattering angle as well as on the energy 
of the photon. 
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Appendices 



A Translation of Tensor Decompositions 

On the left-hand side we give the Passarino-Veltman functions, used in [29l, according to the tensor decomposition 
of Feynman diagrams with respect to external momenta as systematically introduced in lISTl . On the right-hand 
side we follow the corresponding notation in the LoopTools package ll30l . with A 7^ 0. There are also sign 
differences reflecting different notions of metrics. 









(A.i) 






— Ui — O2 


/ A 0\ 

(A.z) 






-L-2 










/A A\ 

(A. 4) 


L'21 




— Oil — ^ 012 — O22 


(A. 5) 






— 022 


(A.o) 


(-^23 




— U12 — O22 


/ A n\ 
(A.7) 






D1+D2+ D3 


(A. 8) 


Dx2 




D2+D3 


(A.9) 


Du 




D3 


(A. 10) 


D21 




Dii + 2 D12 + 2 Z)i3 + 2 D23 + D22 + Ds3 


(A. 11) 


D22 




2D23 + D22 + D3S 


(A. 12) 


D23 




D33 


(A. 13) 


D24 




D12 + D13 + 2 D23 + £'22 + D33 


(A. 14) 


D25 




Di3 + D23 + D33 


(A. 15) 


D26 




D23 + -D33 


(A. 16) 


D27 




-Doo 


(A. 17) 



B Renormalization 

A detailed formulation of renormalization of fermion pair production can be found in various textbooks, e.g. l39l . 
To complete the documentation of our calculation we present some relations resulting from the application of an 
on-mass-shell renormalization, closely following 132 1 . They had been used to derive the formulae given in Section 

HI 

After the renormalization of the boson self-energies: we have to use the following expressions: 

Sr(/) = Sz(/)-ReSz(M|)^Sz(^52)-<5M| (B.l) 
= (B.2) 
E;<f(p2) = s^z(p') (B.3) 

The divergent parts of those renormalized self-energies were given in (13.71) . For the mixing angle renormaliza- 
tion Re i;^(M|) = 5M| and Re ^^/(M^) = are needed: 
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Among the free parameters of the theory we have, only one couphng constant e, using 

g sin 6vi = g' cos = e = 47raem(0) (B.5) 
The electric charge renormalization differs in pure QED and electroweak theory 

(B.6) 

rS^(p') (B.7) 





= 4TTaem{0)ll + 2 


e 




^gweak 

e 


p - 



p2=0 

cos^ M| 



z. 



uv _ e2 1 / 1 1 1 m2,\ i 



sm 



(B.8) 



The wav-function renormalization factor Zf is obtained from the fermion self-energy Ej, with 
The resulting Z factor is: 

Zf = I + ZaJ + ZbJ 15 ■ 

= l + B{m}) + 2mfA'{p^)\^2+C{m})j5. (B.IO) 
For QED, the axial terms vanish, of course. Explicitly, the UV-divergent parts are given by : 



J weak,f/y ;l 

~e " ~6l ^^-^ ^ 

.2aUV _ 2^41 21 2^ 11_4^\ 1 1_2/41 11 2^ M 



<5sin^C' = eU—-—cos'9^ + —cos''9^]-^- = e'{ — -—cos'e^]- (B.14) 
V 6 2 3 / sm^ 9„ e \ 6 3 J e 

with /' denoting the isospin partner of /. 

The above relations define the complete renormalization procedure needed for our reaction. A vertex renor- 
malization, e.g. resulting from terms such as e^j^A^j,"^ in the Lagrangian, traces back to Se and Zf. Explicit 
formulae may be found in the Fortran code (T7i . 

C Infrared Divergences 

The conventions of the one-loop functions and related ones are those used in the package LoopTools l30ll . In 
particular the normalization of the one-loop integration is used as in the following simplest example: 



M^^') = j = -(4V)-^^r(i--) 
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m 



TTl 1 

In ^ + - - 7i? + lii(47r) 
/x^ e 



+ 0{e) 



(C.l) 



In the Fortran program, we leave the treatment of IR divergences to the packages used for the calculation of 
one-loop integrals. Additionally, we checked analytically their cancellation. For this purpose, we isolated them in 
the few IR-divergent scalar integrals contributing to the process e+e^ tt . 

One-loop infrared divergences are due to the exchange of a photon between two massive particles, which occur 
also as external (on-shell) ones. 

Wave function renormalization yields IR-divergent contributions DBq and DBi, the on-mass-shell derivatives 
of -Bo and Bi (w.r.t. the external momentum squared). From ^ With the representation 



1 



one arrives at 



L>Si(m^;m^O) 



2p2 

-^Si(p2;rn2,0)| 



(C.2) 



1 



2m^ 



p =m 

-Ao(m2) +m2 5o(m2;m2,0) - 2171^^0 Bo{p'^;m'^ ,0)\ 



p =m 



(C.3) 



The UV divergences cancel at the right-hand side and the IR divergence is traced back to DBq. We men- 
tion for completeness that the similar function DBi{rTn?\Q,rin?) = -^Bi{p^ = m?]{)^rin?) = [^o("^^) — 
w?BQ{m?; 0, m?)]/ {2m'^), arising from the charged current self-energy with a massless neutrino, is finite. 
An explicit calculation gives 



DBo{p^-rr?M 



1 



r(3 - d/2) 



2 2 



(^2)3-<i/2 (d-3)(d-4) 

Co{m'^, 0, m^; 0, rin?,rin?) 



i2^fi) 



4-d 



d-3 



Co(m^, 0, m^; 0, m^, m^) 



1 m 
— 2 ^'^T- 



(C.4) 
(C.5) 



Assigning the loop-momentum k to the photon line in the initial- and final-state vertex diagrams ensures that 
the divergent part is exclusively contained in one scalar three-point function Cq : 



I 1 2 n 2 2\ 

Co[m ,s,m -,0,171 ,m ) 



s(3 



My) 



1 

21n(l + y) ln{y) - In 



TT 



+ — + 2U2{-y) 



(C.6) 



with 



y = y\s, m, m) = — , ^ = = ^ ^ + ie. 



(C.7) 



Y^l-4m2/s + l + 1 
The finite, small photon mass A is defined according to (I4.70t . 

Finally, IR-divergent functions from the photonic box diagrams, Dq, , D^, and Dfj_^, have to be considered. We 
indicate for the direct two-photon box, shown in Fig. how the singularities can be isolated. 



^In FF, there is no DBi foreseen, while in LoopTools this function was numerically unstable for A 7^ 0. This might be improved 
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+p4 



k-q2 
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Figure .4: An infrared-divergent box diagram. 



The key ingredient for this method 1311 is the following identity: 



Dn oc 



/ 



[{k — — mf] [{k — — mg] [{k — q-^Y — "T-s] [k"^ — w-4 



(C.8) 



2 {k - gi) {k - qs) d^k 



[{k — gi)^ — m-i] [{k — 92)^ — "T-il [{k — 93)^ — [k'^ — ml 



d^k 



[(k - ^2)^ - ml] [{k - qsY - ml] [k"^ - ml] [{k - qi^ - mf] [{k - ^2)^ - mj] [k^ - ml] 

For mi = 7713 = 0, evidently the numerator of the first of the three terms makes it an IR-finite contribution and 
the other two are Cq functions. To demonstrate more exphcitely the procedure we select the above diagram (see 
Fig-Bl and obtain (see also l40l ): 



s Do cx 



IR 



/ 



d'^k 



[{k - qiY] [{k - q2Y - ml] [{k - q^f] [k^ - ml] 
d^k 



+ 



d^k 



IR 



[{k - q2Y - mf] [{k - qs^] [k'^ - ml] J [{k - qiY] [{k - 52)^ - mf\ [k'^ - ml] 

Co{q2, {qs - q2f,qhml,ml, 0) + Co(9i, (^i - ^2)^ 9^; ml, 0, m}) 
2 Co(t, mt,ml;ml,mf, 0). 



(C.9) 



Only one scalar function has to be calculated, and in the limit of vanishing rrie we find 



Co{ml,t,ml;0,ml,mt] 



1 
T 



m^mt . X 
in — — — in 



memt 



+ Lis 1 



iln^ln^ 
2 T T 



(C.IO) 



From the crossed box diagram, we get another function, Dq, with t in Eq. dC.lOl) replaced by u. When combining 
virtual and soft corrections, the singularities of these functions are cancelled against the divergent parts of (14.781) . 
The vector and tensor functions may be treated quite similarly: 



s Dfj, oc 



[{k - gi)2] [(A; - 92)2 _ ^2] p _ ^3)2] [^2 _ ^2 
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IR 



IR 



[{k - - ml] [{k - q^f] [k"^ - ml 



iqi,+qs,)Coimlt,ml,0,mlml), 



+ 



[{k - qiY] [{k - q2Y - rn^Wk"^ - m% 



(C.ll) 



oc 
IR 



/ 



[{k - [(A: - 92)2 _ ^2] p _ ^3)2] [^2 _ ^2] 
{qi^iQiv + q3fiq3u)Co{ml,t,mt;0,ml,mf). 



(C.12) 



To cross-check the result, we isolated the IR-divergent parts also with another approach, where the tensor 
integrals are reduced to scalar ones by means of recurrence relations II411I42II . The divergent contributions hidden 
in the tensor integrals manifest themselves in the form of the three IR-divergent scalar functions Co introduced 
above, namely (IC.5t for self-energies, (IC.6I) for vertices, and dC.lOt for boxes, correspondingly. 
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